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Abstract. 

We consider the propagation of galactic cosmic rays in the 
fractal interstellar medium. Steady state solution of the frac- 
tional diffusion equation, describing cosmic ray propagation, 
is found. We show that the exponent of the steady state spec- 
trum turns out to be equal to the exponent of the cosmic ray 
spectrum above the "knee", that is 77 ss 3.1. 



1 Introduction 

The steady state diffusion equation is frequently used for in- 
terpretation of the cosmic ray phenomena (see, for exam- 
ple,Ginzburg and Syrovatskii (1964); Berezinsky et al. (1990); 
Blumen et al. (1991); Zirakashvili et al. (1991); Ptuskin et 
al. (1993); Kalmykov et al (1999)). Without energy losses 
and nuclear interactions this equation for concentration of the 
cosmic rays with energy E generated by sources distribution 
with density function S(r, E) has the form: 



D(E)AN(r, E) + S(r, E) = 0, 



(1) 



where D(E) is a diffusivity. The equation (1) describes the 
cosmic ray propagation under the assumption that nonho- 
mogeneities of medium have small-scale character. If the 
medium has a fractal structure, the normal diffusion equa- 
tion (1) is not a proper one. A generalization of this equation 
leads to an anomalous diffusion (see the reviews by Bouchaud 
and Georges (1990); Isichenko (1992); West end Deering 
(1994); Uchaikin and Zolotarev (1999) and paper Lagutin 
and Uchaikin (2001)). In recent papers (Lagutin et al., 2000, 
2001; Lagutin and Uchaikin, 2001) the anomalous diffusion 
equations under the different approximations were formu- 
lated and the solutions for point instantaneous and impulse 
sources were found. 

The main goal of this paper is to find the steady state so- 
lution of the fractional diffusion equation and to evaluate the 

Correspondence to: A.A.Lagutin(lagutin@theory.dcn- 
asu.ru) 



exponent of steady state energy spectrum. 



2 Fractional diffusion equation 

The fractional diffusion equation formulated in (Lagutin and 
Uchaikin, 2001) has the form : 

dN{r d f^ = -D(E,a,/3)D^(-Ar^N(r,E,t) 

+ S(r,E,t), (2) 

here D(E, a, 0) = D (a, P)E S is the anomalous diffusivity, 
a, (3 are determined by the fractal structure of the medium 
and by trapping mechanism, correspondingly. Dq + means 
fractional derivative of Rieman-Liouville by time (Samko et 
al., 1987): 



A',' / = T^-^i I (* - T)-»f(T)dT, /* < 1, 



r(i - n) dt 



(— A) Q ' 2 is fractional Laplacian (called "Riss' operator") (Samko 
et al., 1987): 

(_ A W2 f(x) = 1 f A UM d 

R m 

where I > a, x € R m , y e R m , 



d m , l (v)= I {l-e^) l \y\- m - v dy, 

If a = 2 and /3 = 1, the equation (2) is the normal dif- 
fusion equation. If a < 2, /3 = 1, we have from (2) the 



superdiffusion equation: 

^^l = -D(E,a)(-A)^N { r,E,t) 



S(r,E,t), (3) 



discussed in Lagutin et al. (2001). 



3 Steady state solution of superdiffusion equation 

In steady state case the equation (3) takes the form: 

D(-A) a/2 N(r, E) = S(r, E). (4) 

The Green's function G(r, E; Eq) satisfies the equation: 

D(-A) a / 2 G(r, E; E Q ) = S(E - E )6(r). (5) 

The solution of the equation (5) can be found by means of 
Fourier transformation: 

f(k)=Ff(r) = J ' e lkr f{r)dr. 

IP 

Taking into account that the Fourier transform of fractional 
Laplacian is (Samko et al., 1987) 

F{-A) a ' 2 G{r,E-E Q ) = \k\ a G(k,E;E a ), 

it's easy to find: 

5(E - E Q ) 



G(k, E- E 



D\k\ 



S(E-Eo) / dye- Dlk ^ y . 



Applying the inverse Fourier transformation, we obtain: 
G(r, E; E ) = ^-^1 J dy J dfe-^W.. 

_R3 

Since 

9t\r) = ^ f jdke-^-^ 

R 3 

is density of three-dimensional spherically-symmetrical sta- 
ble law, for Green's function we have: 

o 

(6) 

Using Green's function (6) we can find the solution of 
equation (4) for a interesting for astrophysics source. Thus, 
for point source with inverse power spectrum relating to su- 
pernova burst S(r, E) = SoE~ p 6(r), the solution of the 
equation has the form: 

oc 



Taking into account Mellin transform of spherically-symmetrical 
three-dimensional stable laws (Uchaikin and Zolotarev, 1999): 



( a )f a \ - / n( a )(r\ r s - l rl r - 



9^{s)= g£>(r)r'- L dr = 



2*r(j)r(n^) 

a(47r)™/ 2 r (2^i) ' 



we have 



N(r,E) = 



2-*Sp r(^) ^_ 



n 3/2 Dor 3- a r (|) 



£" 



(8) 



(9) 



To evaluate numerically the value of spectral exponent ry = 
p + 6, let's consider in details a passage from non-stationary 
solution to steady state one in a special case a = 1. This 
choice is due to the fact that g^ has the analytical represen- 
tation. It's three-dimensional Cauchy's density: 



9^ 



1 



7T 2 (1 + r 2 ) 2 



Based on the results obtained in Lagutin et al. (2001) for 
points impulse source 

S(r, t, E) = S E- p S(r)Q(T - t)Q(t), 

in our case we have 



N{r,T,E) = 



SqE-p- 5 
2nD r 2 T 



1- 





1 


,,2 


E 2S + 1 



(10) 



Using the representation N(r, E) = NqE v from (10), one 
can easily find the spectral exponent: 



r](T) = p + 6 



26 



TDgE 6 



(11) 



It should be noted that the solution (10) has the "knee" at 

/ \ 1/s 
E (T) = [j%-) ■ At the "knee" energy E the spec- 
tral exponent for observed particles r\ is equal to the spectral 
exponent for particles generated by the sources: 

rj = p. 

One can also see from (11) that at E <C E and E ^> E , we 
have correspondingly: 

VE^E„ ttp-5, 



VE^>Eo ~P + S, 



(12) 



that is the spectrum steepening on 26. The steady state so- 
lution is connected with N(r, T, E) by means of passage to 
the limit: 



So, we have 



N(r,E) = lim N(r,T,E). 

T—*oo 



f]= Km 7](T)=p-\-6. 

T^oo 



(13) 



(14) 
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Fig. 1. Modification of the energy spectrum versus T (Te — > oo). 
Square-position of the "knee" 



E = lim E (T) -> 0. 

T— »co 



(15) 



Thus, the analysis presented above shows that the exponent 
of the steady state spectrum turns out to be equal to the spec- 
tral exponent above the "knee", that is p + S w 3.1. This 
conclusion is illustrated by (fig.l). 

4 Steady state solution of fractional 
diffusion equation(a < 2, f3 < 1) 

It has been shown in (Lagutin and Uchaikin, 2001) that the 
solution of equation (2) for the point impulse source has the 
form: 



N(r,t,E) = 



SnE-P 



D(E,a,l3) 3 / a 



-3/3 /c 



max[0,t-T] 



x¥ 3 a ' 0) (\r\(D(E,a,(3)Tp) 



0\-l/a 



,{a,/3). 



where the scaling function \I/ 3 (r), 

00 



(16) 



(17) 



is determined by three-dimensional spherically-symmetrical 
stable distribution q 3 (r) and one-sided stable distribution 



q[ ' (t) with characteristic exponent (3 < 1 (Uchaikin and 
Zolotarev, 1999): 



7+ioo 



q[^\t) = {2m)- 1 / e^^dX 



(18) 



7— 100 



As the integral diverges if T — > 00, the spectral exponent 
has been evaluated for T ~ 10 10 y. We found r\ w p + 4. 

5 Conclusion 

We have considered the propagation of galactic cosmic rays 
in the fractal interstellar medium. Steady state solution of the 
fractional diffusion equations describing cosmic ray propa- 
gation have been found. We have shown that the exponent of 
the steady state spectrum turns out to be equal to the exponent 
of cosmic ray spectrum above the "knee", that is r\ w 3.1. 
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